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ABSTRACT 

We have carried out a search for galaxies at z ~ 7 - 10 in ^ 14.4 arcmin^ of new NICMOS parallel imaging 
taken in the Great Observatories Origins Deep Survey (GOODS, 5.9 arcmin-), the Cosmic Origins Survey 
(COSMOS, 7.2 arcmin^), and SSA22 (1.3 arcmin^). These images reach 5 a sensitivities of 7no = 26.0-27.5 
(AB), and combined they increase the amount of deep near-infrared data by more than 60% in fields where 
the investment in deep optical data has akeady been made. We find no z > 7 candidates in our survey area, 
consistent with the Bouwens et al. (2008) measurements at z 7 and 9 (over 23 arcmin^), which predict 0.7 
galaxies at z ~ 7 and < 0.03 galaxies at z 9. We estimate that 10-20 % of z > 7 galaxies are missed by 
this survey, due to incompleteness from foreground contamination by faint sources. For the case of luminosity 
evolution, assuming a Schecter parameterization with a typical (p* = 10"^ Mpc~^, we find M* > -20.0 for 
z ^ 7 and M* > -20.7 for z 9 (68% confidence). This suggests that the downward luminosity evolution of 
LBGs continues to z ^ 7, although our result is marginally consistent with the z ^ 6 LF of Bouwens et al. 
(2006, 2007). In addition we present newly-acquired deep MMT/Megacam imaging of the z 9 candidate 
JD2325H-1433, first presented in Henry et al. (2008). The resulting weak but significant detection at /' indicates 
that this galaxy is most likely an interloper at z 2.7. 

Subject headings: galaxms: high-redshift- galaxies: evolution - galaxies: formation 



1. INTRODUCTION 

Populations of Lyman break galaxies (LBGs) have now 
been identified up to z ^ 6, when the universe was less than 1 
Gyr old. Observations now point to earlier times as an impor- 
tant period in the evolution of galaxies. First, some galax- 
ies at z 6 have well established stellar populations, with 
ages ^ 100 Myrs and masses > 10'" Mp (Eyles et al. 2005 ., 
2007; lYan et al.ll2006l: IVerma et al.ll2007riStark et al.ll2007ah . 
requiring significant star formation at z > 7. Second, these 
"first galaxies" likely played an essential role in the reion- 
ization of the intergalactic medium, which occurred some- 
time between z ^ 7 and 1 4 (from the Wilkinso n Microwave 
Anisotropy Probe, WMAP: lOunklev et al.ir2008l) . 

Observations of these z ^ 1 galaxies are crucial; how- 
ever, the search has been significantly more difficult than sur- 
veys for LBGs at z ^ 3-6. At z ~ 7-8, the Lyman break 
passes into the z-band, and galaxies must be identified with 
near-infrared imaging, where sensitivity and area are lim- 
ited. To make matters more challenging, evolution of the UV 
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luminosity function shows declining numbers of luminous 
LBGs with increasing redshift, over the period of z ^ 3 - 6 
jBouwens et al.ll2007i) . Regardless of whether this trend con- 
tinues to z ~ 7, the low density of luminous LBGs at z ^ 6 (a 
few hundred degree"^ to zgso = 26) means that both wide area 
and sensitivity are necessary to continue the search to z > 7. 

Progress in this search for high-z LBGs has been made 
on three fronts. First, wide-area surveys have probed the 
bright end of the luminosity function (LF). Mannucci et al] 
( I2OO7I) used the VLT/ISAAC NIR data in GOODS South to 
search 130 arcmin^ to J ^ 25.5, and Stanwav et all (l2008h 
searched eleven independent sight lines covering 360 arcmin^ 
to Jab = 24-25. Both teams find only a few marginal can- 
didates which they interpret as probable interlopers. Their 
limits are rough ly consistent with e xtrapolation of the z ^ 6 
LF, although .Mannucci et al.l (l2007l) report a slight decline to 
z ^ 7. At higher redshifts, we have searched 135 arcmin^ of 
deep Ji 10 and i/iso parallel images for galaxies a t z 8- 10, 
uncovering one z ^ 9 candidate (JD2325-H1433: iHenry eTaP 
I2008h . 

A second approach has been to use strong gravitational 
lensing to probe the fainter luminosities, where the volume 
density of z > 7 galaxies should be higher Several candi- 
dates have be en found by this technique (Bradley et al. 200j; 
iRichard et a l. 2006, 2 008 ) . How ever, i n an independent anal- 
ysis of the [Richard etsU (l2008h data, iBouwens et"an (l2009h 
suggest that most of these galaxies are either spurious detec- 
tions, or they fail to meet the z > 7 selection criteria. This dis- 
agreement is indicative of the challenge posed by the search 
for these extremely faint galaxies. To make progress, very 
deep observations are needed in both the optical and near- 
infrared. 

This challenge is mitigated by the use of deep NICMOS 
imaging in GOODS, including the Ultra Deep Field (UDF), 
the Hubble Deep Field North (HDFN), and various p arallel 
exposures (.Bouwens et al. 2004b.: 2005, 2008; Labbe et alJ 
l2006HOesch et alJl2009fr Although only eight candidates are 
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Fig. 1. — Cyan squares show the new NICMOS parallel fields in GOODS which we use to search for z > 7 galaxies. Comparison to Figure 1 from Bouwens 
et al. 2008 shows that CDFS-1, -2, and -3 are also used as part of their survey. Coordinates of these fields are listed in Table[T] 
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COSMOS 

Fig. 2. — Cyan squares show the new NICMOS parallel fields in COS- 
MOS, overlaid on the ACS /814 mosaic. For visualization, squares are en- 
larged from the actual NICMOS footprint by a factor of nine in area. While 
fields 4, 5, 8, and 11 are outside the ACS /814 mosaic, they are within the 
Subaru/SuprimeCam z'and i' images. Coordinates of these fields are listed in 
TableO 

found in t hese ^ 23 arcinin^, an d none are spectroscopic ally 
confirmed. lBouwens et al.l (1 200 8') report a luminosity function 
with a bright end that continues to evolve in the same manner 
as those at z ~ 3-6. While uncertain, these data suggest that 
the density of the most luminous z ^ 7 galaxies is even smaller 
than at z ^ 6. 

Because only eight of these galaxies have been found, ex- 
panding the most sensitive combined infrared and optical cov- 
erage to identify even one additional z-dropout LEG would 
be a significant contribution. Accordingly, we have obtained 
~ 14.4 arcmin^ of coordinated NICMOS para llel observations 
in 7iio and Hiau taken in the GOODS fields (Giavalisco et aLl 
2004), COSMOS (Scovilleetal. 2007; Koekemoer et al| 
2007) and SSA22 ( Steidel et al.l 11998). The GOODS and 
COSMOS images reach 5 a = 26.0-26.7 in Juo (0"6 di- 
ameter aperture)- 1-2 magnitudes deeper than the wide area 



gro und based searches ca rried out by iMannucci et alJ (l2007h 
and lStanway et al.l (l2008h . The two parallel fields in SSA22 
are significantly deeper, reaching 5 cr = 27.5 and 27.0 in /no. 
Although most of this area is less sensitive than the UDF and 
HDFN, four out of eight z ~ 7 candidates in Bouwens et al. 
(2008; ^ 23 arcmin^) are bright enough to be detected in 
the deepest of these new GOODS and COSMOS images, and 
most are bright enough to be detected in the SSA22 fields. In 
addition to this search, we have carried out deep follow-up op- 
tical i maging of JD23 25-I-1433, the z 9 candidate presented 
in Hen ry etal.l(l2008h . 

In ^ we describe the data reduction and photometry, as 
well as an overview of the public data products that we use. 
In ^we describe the selection of z > 7 candidates and the cri- 
teria which we use to discriminate against interlopers. In ^ 
we derive a new upper limit on the volume density of z ~ 7 - 8 
galaxies, and discuss implications for the reionization of hy- 
drogen in the intergalactic medium. Finally, in ^we present 
new observations of the z 9 candidate mentioned above, 
which suggest that it is an intermediate redshift interloper 
We use Ho = 70 kms"' Mpc"', = 0.7, = 0.3, and AB 
magnitudes throughout. 

2. DATA 
2.1. Overview 

The data used here consist of NICMOS parallel observa- 
tions taken during GO programs 10872 in GOODS and 1 1236 
in COSMOS (PI H. Teplitz), and 11188 in SSA22 (PI B. 
Siana). For the GOODS fields, 15 parallel fields were ob- 
served in Juo and //leo, and nine lie within the GOODS foot- 
print where ACS data are available. The positions of these 
fields within GOODS are shown in Figure [T] and coordinates 
are listed in Table [T] In total, this corresponds to 5.9 square 
arcminutes of new NICMOS imaging in GOODS. We note 
that three fields (CDF S-1, -2, and -3) are also included in the 
iBouwens et al.l (l2008h search, where they are found not to 
contain any z > 7 candidates. However, in lig ht of the large 
discrep anc y seen in the same NIC MOS data by lRichard et aP 
(2008) and IBouwens et alJ (120091) . we include these fields in 
our search as a consistency check. Typical exposures for these 
NICMOS parallels in GOODS were 8 ks in Jno and 5 ks in 
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The COSMOS parallels consist of twelve fields observed 
in 7no and //i6o', eleven of which lie within the Sub- 
aru/SuprimeCam images in B, r', i', and z' ■ Seven of these 
eleven fields are also within the ACS Isu footprint. A twelfth 
parallel field lies in the north-east corner of COSMOS, where 
the limited SuprimeCam coverage is not sensitive enough to 
discriminate between z > 7 galaxies and interlopers with typ- 
ical galaxy colors. Therefore, we exclude this field from our 
survey. For the remaining eleven COSMOS fields, although 
the optical imaging is not as deep as in GOODS, it is adequate 
to remove interlopers, because, as we will show in § 3.1, no 
z > 7 candidates are found in the COSMOS fields. In total 
these eleven NICMOS parallel fields cover 7.2 arcmin^. Their 
locations are shown in Figure |2] and coordinates are listed in 
Table[T] Typical exposures were 6-8 ks, divided between Jno 
and //i 60. 

Lastly, we include two parallel fields in SSA22, which com- 
prise some of the deepest available NICMOS imaging. How- 
ever, at these faint limits, optical data in SSA22 that are deep 
enough to be useful are limited. Ground-based optical im- 
ages are not sensitive enough to detect the faintest sources in 
the NICMOS images, even if they have typical galaxy SEDs. 
The only available observation that can adequately rule out 
interlopers is an ACS ku image (GO 10405, PI S. Chapman), 
which covers only SSA22-2. Because all NICMOS sources 
are detected in this /814 image, we know that no candidates 
are found in this parallel field (see ® without considering 
z-band data, so we can include it in our survey volume. Al- 
though SSA22-1 can not be used in the z ^ 7 search, we are 
able to use both fields for the /no-dropout LBG search, be- 
cause there are no sources that are red enough in Juo - //leo 
to meet the z ^9 selection criterion in either SSA22 parallel 
field. 

With these data, we select z > 7 candidates as z-dropouts 
and 7iio-dropouts, using the deep optical images to reject in- 
terlopers. This will be discussed in detail in ^ 

2.2. NICMOS Data Reduction 

The NICMOS images were reduced and combined with a 
combination of custom IDL and Python scripts and available 
IRAF procedures. First, images were pedestal-corrected, and 
the South Atlantic Anomaly (SAA) darks were subtracted for 
impacted orbits. Following the SAA correction, the pedestal 
correction was repeated to improve the subtraction. Next, the 
sky frames were made and subtracted using McLeod's NI- 
CRED (1997) code, and a static bad pixel mask which in- 
cluded the vignetted rows was created from these sky frames. 
To remove any remaining gradients in the images, we made 
sky images with each column replaced by its median. This 
image was smoothed by a three-pixel wide boxcar and sub- 
tracted from each NICMOS frame. Then this process was 
repeated for each row of pixels to remove top-to-bottom gra- 
dients. Next, intermittent bad pixels were identified in each 
image using the IRAF package crutil. These masks were com- 
bined wit h the static bad pixel mask, and finally frames were 
drizzled (iFruchter & Hook 1997), using the parameters rec- 
ommended in the dither handbook (pixfrac = 0.6, and scale = 
0.5). Shifts were derived so that the final Ji 10 and i/ieo images 
are drizzled onto the same frame and are therefore aligned. 

' These 12 fields are distinct from the 500 orbits of ffi60 parallel imaging 
in COSMOS (Colbert et al. 2009, in prep), which cannot be used in the z > 7 
search as they lack the essential Jhq imaging. 



The resulting pixels are 0."1, and the zero points that we use 
are adjusted by -0.16 and -0.04 magnitud es in 7i in and i/ isn, 
to correct for the non-linearity reported by'de Jong! ('2006'). 

Sensitivities were measured by randomly placing 0."6 di- 
ameter apertures in the images, rejecting apertures which con- 
tained light from objects^. This procedure is repeated for 
each NICMOS image, as exposure times varied. The 5(j lim- 
its are 26.0-27.5 in Jhq and 25.9-27.0 in //igo, with the faintest 
limits reached in the small area in SSA22 (see Table [U. The 
point-spread-function (PSF) for these NICMOS images was 
derived by stacking several isolated, unsaturated stars. The 
resulting PSF has a FWHM - 0"3 in both bands. The point 
source aperture correction for a 0"6 diameter aperture is 0.31 
magnitudes. 

2.3. Ancillary Optical Data 

GOODS— We use the publicly available v2.0 ACS GOODS 
images in B435, Veoe, '775, and zsso bands-'. Included in v2.0 is 
additional data used to search for Type la supernovae, which 
doubles the vl.O exposure time in the zsso band, and also in- 
creases the sensitivity in (775 . This significantly enhances the 
sensitivity to galaxies at z > 6-7, and improves identification 
of faint interlopers. 

As with the NICMOS images, a PSF is determined by 
stacking several point sources found in the ACS images. We 
find a FWHM of - 0"1 in zsso- Typical 3a Hmits are 28.7, 
28.8, 28.3, in B435, V(,q6, ins, measured in 0."4 diameter aper- 
tures. As we will describe in 32.41 zsso magnitudes are mea- 
sured from 0."6 diameter apertures in images matched to the 
NICMOS resolution. For these, the 3(7 sensitivity is ^27-28 
magnitudes. Some parallel fields near the edge of the GOODS 
footprint have reduced sensitivity. We carefully measure the 
sensitivity in each of the fields, as our objective is to deter- 
mine whether each source is detected in B435, Veoe, or /775. 

COSMOS— The COSMOS data that we use are less ho- 
mogenous than the GOODS data, consisting of both Sub- 
aru/SuprimeCam images at B, r' , /', and z', and where avail- 
able, ACS /814 images. The seeing is 0"8 in B, r', and /', 
and 1"2 in the z' images. Typical limits are 28.4, 27.8, 
27.3, and 26.7 at B, r' , /', and z' in l."2 diameter apertures. 
The ACS /814 images typically reach 27.7 in a 0"4 diame- 
ter aperture. Again, sensitivity varies within the COSMOS 
area, because some of the NICMOS parallels are in the less 
well covered edges. As with the GOODS parallel fields, we 
measure the noise in each field so that we can accurately de- 
termine whether sources are detected in the B, r', /', or /814 
images. 

SSA22 — As described above, the only optical imaging that 
we use for the SSA22 parallels is an ACS hu image that 
covers SSA22-2. We use the "drz" image, directly from the 
archive, which has a 3cr sensitivity of 28.3 in a 0."4 diameter 
aperture. 

2.4. Photometery 

To select z > 7 galaxies, we compare the above described 
z'or Z850 data to NICMOS images. As these data have widely 

^ Apertures containing light from objects were identified in two steps. 
First, we fit a Gaussian to the full distribution of aperture fluxes, including 
those that fell on objects. Then apertures at more than 1 a were rejected and 
the distribution was re-fit. This fit mostly relies on the negative side of the 
flux per aperture distribution. 

' http://archive.stsci.edu/prepds/goods/ 
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Fig. 3. — The two colo r plot allows us to ide ntify z > 7 candidat es from the paralle l s in G OODS (left) and COSMOS (right). The shaded area shows the 
color selection adopted bv lBouwens et al.l 120081) . Model tracks from lBruzual & Charloll j2003l) are: star forming galaxies with E(B-V) =0, 0.5 (solid black, red, 
respectively), a dusty starburst galaxy (red dashed), and an elliptical (black dashed). Numbers indicate fiducial redshift points. 



differing resolution, different techniques are required to mea- 
sure accurate z'- Juo or z%5q -Juq colors. We describe these 
approaches below. 

GOODS — To measure accurate colors of all the galaxies in 
the nine NICMOS fields in GOODs, we downgraded the reso- 
lution of the ACS Z850 images by matching the NICMOS PSF. 
To achieve this, we use the IRAF task PSFMATCH, which 
convolves the ACS images with a kernel made from the NIC- 
MOS and ACS PSFs. The convolved ACS images are then 
rebinned and a ligned with the NI CMOS images. Then, we 
use SExtractor (iBertin & Arnoutsl[T996) in dual-image mode, 
with an inverse- variance weighted average Jwq + H\a) image 
as the detection image. For detection, we require five contigu- 
ous pixels 1 .3 <T above the background. In addition, we use the 
gauss_3.0_5x5.conv filter, which is optimized for finding faint 
sources. Lastly, spurious detections, artifacts, and electronic 
ghosts are manually removed from the catalog. The z^so - Jno 
and JwQ - H\(,() colors are measured in 0"6 diameter apertures, 
and the two-color plot is shown in Figure [3] (left). 

In order to test for non-detections at bands shorter than zs50, 
we measure the flux in 0."4 apertures in the original, uncon- 
volved Z?435, Veoe, and (775 images, at the positions predicted 
by our NICMOS detections. 

COSMOS— The COSMOS data require a different approach, 
because downgrading the resolution of the NICMOS images 
to 1 ."2 seeing causes a significant loss of sensitivity. Instead, 
we resample the z' images to 0"! per pixel (the same as NIC- 
MOS), and ahgn them to match NICMOS. We then used SEx- 
tractor in the same manner as with GOODS, except we use 
0"6 diameter apertures in Jwq and //leo, and l."2 in z' ■ The 
aperture corrections for point sources in these apertures are 
0.3 1 magnitudes for NICMOS and 0.74 magnitudes for COS- 
MOS. Because galaxies at z ^ 7 shou ld be com pact in NIC- 
MOS images dB ouwens et al l l20b4al: iFerguson etani2004t 
iDow-Hygelund et al., 2007h . this treatment is appropriate for 
the sources that we are interested in. For extended sources, 
we expect blue-ward scatter in z' - Juo (away from the z > 7 
selection), as more light will be missed from the higher res- 
olution data. This trend is confirmed for simulated galaxies, 
using the IRAF artdata package. The two-color plot for the 
COSMOS fields is shown in Figure |3] (right). 



As with the GOODS data, we measure the flux at the po- 
sitions predicted by the NICMOS detections, using 0"8 aper- 
tures in B, r' and and 0."4 apertures in /814. 

SSA22 — As we are not using any z- band data for SSA22, 
there is no need to properly account for z- J\ 10 colors mea- 
sured with mismatched apertures and resolutions. Therefore, 
we simply follow the same procedures described above for the 
GOODS and COSMOS parallels- measuring Ji 10 and //160 
with SExtractor, and testing for /814 detections in 0."4 diame- 
ters apertures in SSA22-2 where the ACS data are available. 

3. RESULTS 
3.1. Selection of z > 7 candidates 

Candidates for z > 7 galaxies are selected using the follow- 
ing criteria: First, we require that galaxies are detected in the 
7110 + 7^160 detection image at > 5a significance. In total, we 
find 696 sources that meet this criterion in GOODS, 701 in 
COSMOS, and 21 1 in SSA22. Next, z > 7 candidates must be 
undetected at the 2 a level in bands bluer than z' or Z850- This 
eliminates the vast majority of sources, with only two candi- 
dates remaining in the GOODS fields, two in COSMOS, and 
none in SSA22-2. We list these sources in Table |2] and we 
will proceed to show that all are interlopers. 

We next use the colors of these sources to determine if any 
have SEDs consistent with z 7 galaxies. We adopt the color 
cuf* of Bouwens et al. (2008, 2009), so that candidates must 
have z- /no > 0.8 andz - Jim > 0.8 + 0.4 (Jiur ^^leo), and /no 
-//160 < 1.2 (where z refers to both zgso and z')- All four of 
the "dropout" sources mentioned above lie outside this selec- 
tion. One source (C5-zDl), has z' - /no ~ 0.2, and the oth- 
ers (CDFS3- JD1, CDFS4 -JDI and C8-JD1) have Jno - //160 
> 1 .2. While lOesch et al.l ('2009) have suggested a stricter cut 
of Z850- Ji\o> 1-3, adopting this cut would make no difference 
in our search, because we have not found any candidates with 
the most generous selection. 

The red Juo - H16O colors of these three sources could be an 
indication of the Lyman break in the 7no band, and redshifts 
z > 8 ( Bouwens et al. 2005 ; Henry et ali,2007. 2008). How- 
ever, 7no - dropouts must also be undetected at the 2a level in 

Despite differing filter set for the COSMOS data, this color cut selects 
galaxies at z > 6.5 in both cases. We will show in |4]that the survey volume 
is not affected by this inhomogeneity. 
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the z' or zg,5o bands. This restriction eliminates CDFS3-JD1 
and CDFS4-JD1. The remaining source, C8-JD1, can not be 
ruled out on the basis of optical/NIR data alone, but longer 
wavelength data from IRAC on the Spitzer Space Telescope 
show strong detections at 3.6 and A.Sfim (//leo - [3.6] = 2.4, 
//i6o - [4.5] = 2.7). For z ^ 9, these colors correspond to a 
rest-frame UV slope which is much redder than Lyman break 
galaxies, so this galaxy is more likely an interloper at z 1-3 
with a dusty starburst or an old stellar population. In conclu- 
sion, none of the four optical "dropout" sources that we find 
can be described as a plausible z > 7 galaxy. 

3.2. On Incompleteness from Foreground Contamination 

We have rejected as interlopers any sources which have 2a 
detections in ^435, B, V(,o6, r', i' or /814. While this approach 
is commonly taken in LBG surveys, it does not consider the 
possibility that a weak detection in any or all of these "veto" 
bands could arise from foreground contamination. In fact, as 
we will demonstrate, the probability of contamination is sig- 
nificant. 

To estimate the influence of foreground contamination, we 
use the UDF ACS catalogs. The surface density of sources 
brighter than our typical 2a detection limits in GOODS (B435, 
Vm and /775~ 29.1, 29.2, and 28.7) is - 400 arcmin'l This 
corresponds to about a 10% probability of a foreground con- 
taminant lying within 0."5 of a NICMOS detected source. For 
the COSMOS fields, the B, r', and /' limits ai-e shallower, so 
the surface density of possible contaminants is lower 250 
arcmin"^). However, the seeing-limited resolution requires 
larger apertures. In this case, we find that the probability of 
a foreground source lying within 1" of a z > 7 candidate is 
about 20%. We therefore estimate that 10 - 20 % of true z > 7 
galaxies would be rejected by our survey because of faint fore- 
ground contaminants. 

4. DISCUSSION 
4.1. The z>l Luminosity Function 

While we have not detected any candidate z > 7 galaxies, 
we can place limits on the luminosity functions (LFs) of z- 
dropouts at z ~ 7, and J\u) -dropouts at z ~ 9. Furthermore, 
we can compare this limit to predictions from LFs at z 3-6 
and place constraints on evolution from z ^ 6 to 7. 

First, we calc ulate the survey volume following 
ISteidelet al.l([T999h : 

Veff(M) = JpiM,z)^dz. (1) 

The quantity p(M,z) is the probability of both detecting a 
source of a given absolute magnitude and redshift, and se- 
lecting it as a z- or Jhq- dropout based on the criteria that 
we established in |3] This probability can be expressed as 
p(M,z) = S(M,z) X C{m), with S(M,z) representing the selec- 
tion function, and C(m) the photometric completeness. We 
use simulations to determine these quantities for both z ~ 7 
and z ~ 9. To obtain S{M,z), we require only one assumption, 
namely, a distribution of galaxy spectra. We use a Gaussian 
distribution of UV power-law slopes estimated from z ^ 6 
galaxies (/a cx ; /? = -2.2 ±0.2; Stanway et al. 2005). 
Then, for every M, z, and /3, we predict the z', Z850> Jiio and 
HidQ- band magnitudes, as well as the the magnitude in the 
/no + //16O image that we used for detection. Sources are 
required to be (1) bright enough to be detected at > 5a sig- 
nificance in the Jhq + Hk^ image, and (2) meet the color se- 
lection criteria discussed in ^ We also assume that 15% 




Fig. 4. — The effective volume as a function of absolute UV magnitude 
for the 14.4 arcmin^ covered by this survey (black), and the 5.8 arcmin- of 
the NICMOS UDF (red). The solid curves are for z-band dropouts at z ~ 7, 
and the dashed curves are for 7iio- dropouts at z ~ 9. SSA22-1 is excluded 
from the 7, ^1 search, as it does not have adequate optical data to rule out 
interlopers. The UDF volumes are estimated in the same way as the volume of 
this survey, assuming a photometric completeness similar to SSA22-1, which 
has the same sensitivity as the UDF. 

of all z > 7 galaxies are missed because of foreground con- 
tamination, as we showed in ^3.21 We calculate S{M,z) for 
both the z' (COSMOS) and zsso (GOODS) filter sets, and find 
that the difference is less than 2% (for a fixed 7no + Hk^ 
apparent magnitude limit). Therefore, the only difference in 
the GOODS and COSMOS portions of this survey is that the 
NICMOS images in COSMOS are slightly shallower. 

We measure the photometric completeness, C(m), using the 
IRAF package, artdata, to add point sources to the /no + i/ieo 
images. We then use SExtractor with the sam e con figuration 
that we used for the photometry described in 32.41 We find a 
typical completeness of 70-80% at the 5ct detection threshold 
for the aggressive SExtractor parameters that we have cho- 
sen. Finally, to evaluate Equation (1), we assume that all of 
the z-dropouts are at z = 7, and the J[ 10-dropouts are at z = 9, 
so that C{m) translates to C{M). The resulting effective survey 
volumes for z ~ 7 and z 9 are shown in Figure |4] As men- 
tioned in ^ due to limited optical data, we can only include 
SSA22-2 in the z 7 search, but both SSA22 fields are in- 
cluded in the z 9 upper limit, as they contain no candidates 
with7iio-//i6o > 1.2. 

We next constrain the UV luminosity function. Assuming 
a Schecter parameterization of the LF, we show the space al- 
lowed for (f)* and M* in Figure |5] The shaded areas show the 
upper limits for 68 and 95% confidence for the z ~ 7 survey, 
and the area below and to the right of the dotted lines indi- 
cate the same for the z 9 search^. We also plot measured 
LFs from Bouwens et al. (2007, 2008) at z 4,5,6,7, and the 
upper limit at z ~ 9. The non-detections that we find in this 
survey are consistent with the Bouwens et al. measurements, 
which predict 0.7 z 7 candidates in our survey, although 
the error bars on their z 7 LF are large, due to the small 
sample. The dashed line shows the upper limit at z 7 from 
Mannucci et al. (2007). Their constraint on luminous M* is 
stronger than what we have measured here, due to their wide 

Uncertainties given here are in the Poisson noise limit, which is the 
dominant source of uncertainty when the expected density of sources is < 1 
arcmm l iTrenti & Stiavellill2508h . Cosmic variance is also greatly reduced 
because of the large number of independent sight lines that we have searched. 
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Fig. 5. — Our non-detection of z > 7 galaxies constrains the luminosity 
function of these galaxies. The shaded areas indicate the allowed area for 
M^fy and 0* , for 68% (light grey) and 95% (dark grey) confidence. Atz^9, 
these same upper limits are shown by the dotted lines, with the allowed pa- 
rameter space down and to the right. Here, we assume a Schecter parameter- 
ization of the LF, with a faint end slope of a =-l .74 ( Bouwens et al. 2007). 
The redshift labels refer to the measurements from Bouwens et al. (2007, 
2008), marked by the points, and the aiTow indicates their M* upper limit 
at z ~ 9 when (p* = 10"' Mpc"'. The square is the LF measurement from 
lOesch et al.l <2059l) , which holds (/>* fixe d at 1.4xlO~' ^P^'!: '^^^ dashed 
line indicates the upper limit (68%) from lMannucci et al] j2007h . again, with 
the allowed parameter space down and to the right. 

area survey (^130 arcmi n^). Our result i s also consistent with 
the constraint reported by Stanwav et al.l (12008 ). where a z ~ 7 
upper limit that is similar to the z ~ 6 LF is found. 

This limit can be used to address the controvers y over 
the numbers of strongly len sed z > 7 galaxies (Richa rd et al.l 
I2008t iBouwens et al.ll2009h . These authors have found dif- 
fering numbers of candidates behind the same lensing clus- 
ters, using the same NICMOS data. Richard et al. find a few 
times more candidates than are predicted from the small un- 
lensed sample in the field (Bouwens et al. 2008). In fact, such 
a comparison is difficult to make, as the lensed and field sur- 
veys observe mostly different ranges of luminosities. While 
the unlensed < J+H > apparent magnitudes of the Richard et 
al. sources range from 27-30, the Bouwens et al. field survey 
finds sources down to //if,o ^ 28. However, within this one 
magnitude of overlap, the Richard et al. density agrees more 
with the Bouwens et al. measurement at z ^ 6 than at z 7. 
While our survey probes even brighter magnitudes, we can 
compare to the Bouwens et al. 6 LF. Assuming no evolu- 
tion, this LF predicts 3.2 z ~ 7 galaxies in our survey volume- 
a scenario which we can exclude with 97% confidence (Pois- 
son statisti cs). Our result is mor e consistent with the z ^ 7 
result from lBouwens et aTl ( l2008h . as shown in Figure|5] 

4.2. Star Formation and Reionization 

We also constrain the amount of star formation at z ^ 7 
and 9. To do this, we fix </)* ~ 10"^ Mpc"^ mag"'. This 
is supported by lumino sity functions that have been mea- 
sured by many authors (IBouwens et al.l l2007, and references 
therein), from z ~ 3 - 6. While some scatter is present at 
z ~ 6, most LFs agree with this value of (j)* to within a fac- 
tor of two, so that any evolution of this parameter must be 
small. For this choice of (j)*, we find that M* > -20.0 at 
z ~ 7, and > -20.7 at z ^ 9. Assuming a steep faint end 
slope of a = -1.74 (IBouwens et al.l 120071) . and integrating 
the LF to zero luminosity, we find a luminosity density of 



PL < L5 X 10-'' ergs-'Hz-'Mpc"^ at z 7. This limit is 1.9 
times higher at z ~ 9. This corresponds to a star formation 
density of pspR £ 0.019 Mfriv r "' Mp c"^' at z ~ 7, when the 
conversion from iMadau et al.l (1 19981) is used. It is impor- 
tant to note that this conversion assumes no extinction, solar 
metallicity, and a Salpeter IMF with dN/dM oc M"^^ from 
M = . 1 - 1 00 Mq . While a correction to a more likely metal- 
licity of 0.2Zo is negligible (< 5 %), a shallower IMF slope 
of -1.7 will decre ase the SFR by a factor of 3.2 (calculated 
from Starburst99: lLeitherer et anil999l) . 

An important question remains whether galaxies at z 6-7 
are capable of reionizing the neutral hydrogen in the inter- 
galactic medium. This question is difficult to address, as it 
depends on the duration of the reionization. A longer reion- 
ization will require more ionizing photons over the li fetime 
of the galaxies in order to account for recombination ( Char yl 
|2008|) . Nonetheless, it is interesting to compare our upper 
limit to the recombination rate at z 7, for a completely ion- 
ized IGM (consist ent with the WMAP 5 year electron s catter- 
ing optical depth; iDunkley et ani2008h . iMadau etalJ (1 19991) 
report this rate in terms of the critical SFR required to main- 
tain an ionized IGM: 

_ 0.039 Mq yr-i Mpc'^ / 1 -i-z\ V C \ / ^hh^ \ ^ 

(2) 

where, again, solar metallicity and a Salpeter IMF from 0.1- 
100 Mq are assumed. This critical SFR also depends on a 
number of other important, but uncertain parameters. The 
escape fraction of ionizing photons, fesc, has been difficult 
to measure. While a number of authors have found that 
the escape fraction is small (< 5-10% relative to photons 
escaping at 1500 A'': iMafkan et al.1 llOOl ISiana et al.ll2007l 
Bridge et al., in prep), ther e remains some evidence that it 
could i ncrease with redsh ift (ISteidel et al.l2001l: IShaplev et al.l 
I2006t llwata et al.1 [2008). The HII clumping factor, C = 
("ifn) / i'^Hii)^ , is also important, as this dictates the average 
recombination rate per hydrogen atom relative to an IGM of 
uniform density. While many authors have adop ted an esti- 
mate o f C = 30, based on simulations by Gnedin ^ OstrikeJ 
(Il997l) . more recent work suggests that this estimate is 
much too high, and C < 10 may be more appropriate (e.g. 
iBolton & Haehne lt 2007; Trac & Cenll2007h . 

In order to meet the requirement posed by our upper limit 
of 0.019 Moyr"' Mpc"^ at z 7, we find that C/fesc < 15. 
However, this number is strongly influenced by the faint end 
slope of the LF, because we have integrated our constraints 
to zero luminosity. We have assumed a f aint end slope of 
a = - 1.74, based on IBouwens et alj (|2007|) . but lOesch et alj 
(l2007l) show that this slope is influenced by input assumptions 
such as dust extinction and IGM neutral hydrogen absorption 
(which alter the effective survey volume). For the shallower 
slope of a = -1.6, reported by Oesch et al., our upper limit 
is reduced by a factor of 1.7, and we then require C/ fesc ^ 9. 
On the other hand, it has been predicted that a approaches -2 
for a sample of young galaxies undergoing their first signif- 
icant bursts of star formation (Overzier et al. 2008). In this 
case, constraints are more dependent on the true low luminos- 
ity cutoff. 

* These escape fraction upper hmits from the literature are the relative es- 
cape fraction, described by Shapley et al. (2006) and Siana et al. (2007), as 
opposed to the absolute escape fraction that we use in this paper. By defini- 
tion, the absolute escape fraction is smaller than the relative escape fraction. 




J^PJ^pW ;f160W 3.6 |xm 4.5 ^im 

Fig. 6. — Postage stamp images of JD2325+1433, from left to right: /'(Megacam), Jno, //i60 (NICMOS), 3.6fim, and4.5/im (IRAC). Images are 7 "on a side, 
and are oriented witli nortli up, and east to the left. Exposure times in Jhq and 3.6/jm are several times longer than their //i6o and 4.5/im counterparts, so the 
photometry does indeed suggest two spectral breaks. The NICMOS and IRAC images are described in more detail in Henry et al. (2008). 





The effects of metallicity and IMF are also important in de- 
termining the ionizing output of galaxies. We use Starburst99 
models ( Leitherer et al. 1999) to calculate the ionizing photon 
rate for metal poor stellar populations (Z = Q.2Zq) and for a 
shallower IMF slope. For a Salpeter IMF and Z = Q.2Zq, a 
stellar population will produce 1 .4 times more ionizing pho- 
tons than a solar metallicity population with the same UV lu- 
minosity. Consequently, the constraint from this survey be- 
comes C/fesc < 21. Likewise, with Z= O.2Z0 and a shallower 
IMF slope of dN/dM oc M~'^, this constraint is relaxed to 

C/fesc < 36. 

Lastly, it has also been noted that the electron temperature 
in the primordial HII regions will play an i mortant role (e.g. 
iTumhnson et al.ll2001l: IStiaveni et 311120041) . Because the re- 
combination coefficient is proportional to 7""°^, a factor of 
two increase in temperature decreases the critical star forma- 
tion rate by a factor of 1 .6. 

In summary, we find that for reasonable models, C/fesc 
30-40 is required to maintai n an ionized IGM at z ~ 7. This 
echos constraints reported by , Chary! (l2008h . who finds that for 
C/fesc ~ 60 ("high-V" case) and a Salpeter IMF the number 
of ionizing photons produced is too low to reionize hydrogen, 
unless the reionization occurred rapidly between 6 < z < 7. 

5. FOLLOWUP OF THE Z ~ 9 CANDIDATE JD2325+1433 



In iHenrv et"al] (l2008l) we reported the discovery of a 
luminous z ~ 9 candidate from the wide area, NICMOS 
Pure Parallel Su rvey ( 135 arcmin^ to 7|io and //ifip ~ 25 
AB: iTephtz et al.l lT99l lYan etal.ll2000t IColbert et al.ll200l 
iHenrv et al.ll2007l) . This candidate, JD2325H-1433, was iden- 
tified as having a strong spectral break between the Ji lo and 
i/i6() bands, with a faint but detected /no flux and /no - //leo 
= 1.7. Subsequent followup observations with SpitzerflRAC 
showed a flat spectrum in //leo - [3.6], and a second spectral 
break between 3.6 /im and 4.5/im. The only possibility for 
two breaks are the Lyman and Balmer breaks, and a redshift 
of z ~ 9. However, given the uncertainties in IRAC flux, the 
significance of the second break is only about 95%, and with- 
out this break, the galaxy spectrum could also be fit by an 
intermediate-redshift elliptical or post-starburst galaxy. 

The main impediment to a robust identification of 
JD2325-I-1433 as a z 9 galaxy is the lack of deep optical 
imaging to verify that we have indeed identified the Lyman 
break. Such observations require a significant investment, and 
a non-detection at / 28 AB would ultimately not be defini- 
tive because interlopers could be even fainter than this. On the 
other hand, obtaining a detection would definitively rule out 
the z 9 interpretation. Therefore, we have obtained /' obser- 
vations with the MMT to attempt to understand the nature of 
JD2325H-1433. 

5.1. MegaCam Observations of JD2325+1433 



The i' observations consisted of a series of exposures, of 
length 300 to 500 s each 6.8 hours), taken on the nights of 
2008 June 19-24 with Megacam at the 6.5 m MMT (Mcleod 
et al. 2006). The observations were carried out through thin 
cirrus, except for the nights of 2008 June 20 and 21, which 
were photometric. Seeing varied from as low as 0"8 to as 
high as 1 "6, and averaged about 1 ."0. The data were reduced 
interactively using standard techniques: bias-subtracted and 
flattened exposures were treated to remove cosmic rays and 
bad pixels before calculating the coordinates using stars from 
the USNOBl.O catalog and correcting the photometry for off- 
axis scattered light. The resulting exposures were spatially 
registered to a common coordinate system. All frames taken 
on the nights of 2008 June 19, and 22-24 were then flux- 
calibrated using exposures from the photometric nights, and 
all the frames were then coadded to create an /' mosaic. The 
final image has a seeing FWHM of 1"2. We show a cutout 
image centered on JD2325-I-1433 in Figure |6] a longside our 
NICMOS and IRAC images that are described in lHenrv et alj 
(.2008.) . 

5.2. Photometry 

We use the NICMOS images to predict the position of 
JD2325-I-1433 in the i' image, and measure the flux in a l."3 
diameter aperture at this position. The noise is measured by 
randomly placing apertures in blank parts of the image, as was 
done with the NICMOS and other optical images (see We 
find a S/N of 2.6, and an aperture magnitude of 26.8 ± 0.4. 
The aperture correction measured for point sources in the field 
is 2.18 ±0.04 in flux units, and so the result is /'= 26.0 ±0.4, 
total. Although the detection is weak, it strongly suggests an 
intermediate redshift interloper. The probability of the /' de- 
tection being the result of a foreground contaminant (as we 
described in within 1" of JD2325H-1433 is low (- 5%), 
as the /' image is not as deep as the GOODS and COSMOS 
optical images. 

5.3. An Updated Photo-z of JD2325+1433 

We update the photometric redshift of JD2325-I-1433 by in- 
cluding th e i' measurement, a nd repeating the fit that we per- 
formed in iHenrv et"an (1200 8h . To do this, we use the pho- 
tometric redshift code, H yperz (iBolzonella et al.ll2000l) . with 
iBruzual & CharlotI (l2003l) stellar synthesis templates. We fit 
for redshift, allowing age, extinction (using Calzetti et al. 
2000), and metallicity to be free (Z = 0.02,0.2,0.4, and 1 x 
Zq), and using four star formation histories: an instantaneous 
burst, a constant SFR, and two exponentially declining star 
formation histori es with e-fol ding times (tsfr) of 100 and 
500 Myrs. As in iHenrv et all (2008), we do not include the 
upper limits at 5.8 and 8.0 /im, as they do not constrain the 
fit. The revised, best-fitting model is shown in Figure |7] It is 
described by a 250 Myr instantaneous burst at z ^ 2.7, with 
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Fig. 7. — i-e/f- The addition of /' improves our pliotometric redshift, and z ^ 2.7 is now favored. The solid line is the prefer r ed fit, wiiich is a 250 Myr old 
instantaneous burst model, with Ay = 0.2 and solar metallicity. The grey dashed line is the best fit z ~ 9 SED from lHenrv et al.l 420081) , which is a 64 Myr old 
constant star forming model with Ay = l.O and Z = 0.2Zq . Right- The redshift probability distribution from our Monte Carlo simulation shows a peak at z ~ 3 
when the newly acquired /' detection is included. 



solar metallicity, Av/ = 0.2, and a stellar mass of 9.9 x 10^ M0. 
The absolute B-band magnitude is Mb = -21 .0. 

We use Monte Carlo simulations to assess this undercon- 
strained problem by constructing a five dimensional (z, age, 
Av, metallicity, and star formation history) probability den- 
sity function. This is done by generating 10^ realizations 
of the photometry, with magnitudes simultaneously perturbed 
according to the uncertainties. We then repeat the fit described 
above. The probability distribution in redshift space is shown 
in Figure |2l Now, z ^ 2-3 solutions are favored, with 74% 
of realizations having a best fit at z < 5. The fact that the 
Z ~ 8 - 10 interpretation still comprises a significant fraction 
of the simulated fits is guaranteed by the low S/N /' detection, 
which frequently dips below 1 cr when perturbed in the Monte 
Carlo simulation. For these cases we do not include the /' 
observations and the fits strongly favor the z^9 solution. Re- 
gardless, the inclusion of this weak detection in our analysis 
adjusts the preferred redshift to z ^ 2 - 3. This is more in line 
with an extrapolation of the Bouwens et al. (2006, 2007) lu- 
minosity functions, which imply a low likelihood of a galaxy 
at z 9. 

The additional constraints from our Monte Carlo simulation 
suggest, for z < 5: (1) a poorly constrained age with a median 
of 360 Myrs, and a 68% confidence interval ranging from 100 
Myrs to 1 Gyr and (2) little or no extinction, with 68% of 
realizations preferring Ay of 0.5 or less. 

5.4. Interlopers in future z ^ 9 surveys 

The discovery that JD2325H-1433 is an interloper has im- 
portant implications for future z ^ 9 surveys, because similar 
sources will be readily discovered with new near-infrared in- 
struments. In addition to JD2325H-1433 in the NICMOS Pure 
Parallel Survey we find 12 more galaxies down to H\a) ^ 24 
which have similarly red /iio - //i6o> 1.7. As this wide area 
survey is complete for such red galaxies at this limit, the den- 
sity of these objects is approximately 200 degree"^. Longer 
wavelength IRAC observations of a few of these sources indi- 
cate rising SEDs that are indicative of interlopers, but not all 
of these unusually red galaxies have yet been observed with 
IRAC. So it is likely that more galaxies with extremely red 
J no - H\a) and a flat spectrum at longer wavelengths have 
been detected in the NICMOS pure parallel imaging. These 



sources will be difficult, if not impossible to distinguish from 
z > 7-8 galaxies in future surveys, meaning that deep optical 
imaging or a high S/N detection of the B aimer break will be 
crucial. 

Using deep optical imaging to distinguish z > 7-8 galaxies 
from interlopers will be challenging. At the faint magnitudes 
where these galaxies are more likely to be confirmed {H > 28 
AB), optical observations with the James Webb Space Tele- 
scope will take at least 10 hours per pointing to reach > 30 
AB at the 2 a level required for non-detection. In total, this 
investment in telescope time simply to confirm non-detections 
could amount to hundreds of hours. In addition, as we showed 
in ^3.21 foreground contamination from faint, lower-redshift 
objects can be a substantial source of incompleteness. Extrap- 
olating number counts from the UDF, we estimate the surface 
density of galaxies brighter than 30 AB (total magnitudes) 
in B435, Veoe, and i-m is ^ 900 arcmin"^, or 0.25 arcsec"^. 
Clearly, high angular resolution will be necessary to distin- 
guish interlopers from z > 8 galaxies, as ground-based seeing 
limited observations would suffer from severe confusion. 

6. CONCLUSIONS 

The absence of any z > 7 galaxies in our new NICMOS 
data strongly constrains the volume density of z > 7 galax- 
ies. We have shown that at z 7, if (f>* = 10"^ Mpc~^^, then 
M^y > -20.0, and the cosmic star formation density (inte- 
grated to zero luminosity) is < 0.019 Mq yr"' Mpc"^. Al- 
though the luminosities that we observe are much brighter 
than the candidates reported from lensing surveys (Richard 
et al. 2006, 2008), we can indirectly addr ess their discrep- 
ancy with the field survey of lBouwens et alj ( I2OO81) . Our non- 
detection is consistent with Bouwens et al., so our indepen- 
dent result supports their reported evolution for the most lu- 
minous sources. This suggests an additional fading of M^y 
by 0.4 magnitudes at from z ^ 6 to z 7. 

Clearly, large uncertainties remain as the few reported can- 
didates are hardly robust detections. Upcoming surveys us- 
ing the Wide Field Camera 3 (WFC3) on boai'd HST will ad- 
dress this issue with its improved resolution and sensitivity, 
increasing the number of known z 7 candidates by an order 
of magnitude. Current plans to use pure parallel mode obser- 
vations to cover a wide area (Pis M. Trenti, H. Yan, and M. 
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Malkan) will also provide crucial measurements of the lumi- 
nous sources. 

Interpretation of the UV luminosity function in terms of the 
ionizing photon budget required for neutral hydrogen reion- 
ization is uncert ain, for reasons that we (in ® and many oth- 
ers (e.g. Bunker et al.l 12004: Bou wens et al. 2007) have dis- 
cussed. However, for a Salpete r IMF and a faint end slope of 
a = -1 .74 (reported at z 6 bv lBouwens et al.ll2007h . we find 
C/fesc < 15 is required to maintain a com pletely ionized IGM 
at Z^ l ■ Fo r current estimate of C ^ 10 dBolton & Haehnell 
l2007y Trac & Cen 2007) and the commonly adopted fesc = 0.1 
(e.g Ichary 2008|)^ this ratio is C/fesc = 100- far too high for 
star forming galaxies to maintain a completely ionized IGM 
at z ~ 7. However, what is more likely is that our result pro- 
vides indirect evidence for significant evolution in one or both 
of C and fesc- 

We also pres ent foUowup observ ations of the z^9 candi- 
date reported in^ Henrv et al.l ( l2008l) . With deep imaging from 
the MMT we find a 2.6(t detection at which suggests an 
intermediate-redshift interloper. This interpretation of the for- 



mer z^9 candididate, JD2325H-1433, is more consistent with 
upper limits reported by Bouwens et al. (2005, 2008,2009), 
as well as the upper limit which we find in this study. The 
fact that this interloper has such an extremely red J\\q - Hia) 
and LBG-like SED at longer wavelengths means that similar 
sources at fainter magnitudes will require a large investment 
in optical imaging in future surveys, such as those with WFC3 
and in the longer term, JWST and future thirty-meter class 
telescopes. 
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TABLE 1 
GOODS Fields 



ID 


RA (J2000) 


Dec(J2000) 


■6435'' 


^606" 


'775" 4l4 


2850*' 


■fiw'' 




CDFS-1 


03 32 26.78 


-27 41 58.4 


28.7 


28.9 


28.3 


28.0 


26.7 


26.4 


CDFS-2 


03 32 40.08 


-27 44 03.1 


28.6 


28.8 


28.3 


27.8 


26.7 


26.4 


CDFS-3 


03 32 24.31 


-27 40 23.3 


28.7 


28.9 


28.3 


28.0 


26.7 


26.4 


CDFS-4 


03 32 52.30 


-27 46 50.9 


28.7 


28.5 


28.2 


27.7 


26.5 


26.3 


CDFS-5 


03 32 27.92 


-27 40 14.8 


28.8 


28.8 


28.2 


27.9 


26.5 


25.9 


HDFN-1 


12 37 24.56 


62 16 23.9 


28.6 


28.8 


28.4 


28.0 


26.5 


26.1 


HDFN-2 


12 36 06.10 


62 12 16.1 


28.7 


28.9 


28.4 


28.0 


26.4 


26.2 


HDFN-3 


12 27 11.00 


62 15 57.9 


28.6 


28.9 


28.4 


28.0 


26.5 


26.4 


HDFN-4 


12 36 09.09 


62 06 34.1 




28.5 


27.9 • • • 


27.7 


26.5 


26.4 



COSMOS FiclcK 



ID 


RA (J2()()()) 


Dec (J2()()()) 


B" 




,/d 


/814" 




J\ 111' 


W 1 (,()^ 


COSMOS- 1 


10 01 58.51 


02 09 35.6 


28.5 


27.8 


27.3 


27.8 


26.7 


26.4 


26.0 


COSMOS-2 


10 00 32.63 


01 59 23.0 


28.8 


28.1 


27.7 


27.3 


26.9 


26.3 


26.2 


COSMOS-3 


10 00 30.00 


02 02 00.0 


28.8 


28.1 


27.7 


27.8 


26.8 


26.0 


26.2 


COSMOS-4 


10 01 47.25 


02 56 48.2 


28.0 


26.9 


26.7 




26.0 


26.1 


26.2 


COSMOS-5 


10 02 24.91 


02 51 46.3 


28.5 


27.8 


27.3 




26.7 


26.0 


26.1 


COSMOS-6 


09 58 32.62 


01 48 24.0 


28.5 


27.8 


27.2 


27.5 


26.6 


26.2 


26.3 


COSMOS-7 


09 58 27.51 


02 18 29.1 


28.4 


27.8 


27.3 


26.6 


26.8 


26.5 


26.4 


COSMOS-8 


09 58 40.25 


02 52 52.7 


28.4 


27.8 


27.3 




26.6 


26.1 


26.0 


COSMOS-9 


10 02 10.58 


01 45 46.1 


28.4 


27.9 


27.4 


27.7 


26.7 


26.5 


26.0 


COSMOS-10 


09 58 22.61 


02 39 02.3 


28.4 


27.9 


27.3 


27.7 


26.6 


26.1 


26.3 


COSMOS- 11 


10 00 02.82 


02 46 07.8 


28.4 


27.8 


27.3 




26.6 


26.0 


26.2 



SSA 22 Fields 



ID RA (J2000) Dec (J2000) W Jno" Hm" 



SSA22-1 22 17 21.23 00 24 09.8 ■■• 27.4 27.0 

SSA22-2 22 17 23.36 00 22 03.6 28.3 27.0 26.5 



Note. — Sensitivities were measured by randomly placing apertures in blank parts of the images. 
All limits are in aperture magnitudes, and aperture corrections are given in §2. 

^ 3 cr limits measured in 0"4 diameter apertures.'' 3 a limits measured in images that were PSF- 
convolved to match the NICMOS resolution, using ()"6 diameter apertures. 5 cr limits measured in 
0"6 diameter apertures.'' 3 a limits measured in 0"8 diameter apertures."^ 3 a limits measured in 1"2 
diameter apertures. 



TABLE 2 
Optical Dropout Sources 



ID 


RA (J2000) 


Dec (J2000) 


z-Hi(,o 


■'110 - ff|60 




CDFS3-JD1 


03 32 23.24 


-27 40 20.8 


1.7 


1.6 


25.1 


CDFS4-JD1 


03 32 51.66 


-27 47 15.3 


1.7 


1.5 


25.3 


C5-ZD1 


10 02 24.47 


02 52 05.4 


0.2 


-0.1 


25.7 


C8-JD1 


09 58 39.07 


02 52 53.6 


> 1.2 


> 1.6 


25.1 



Note. — //iso magnitudes are aperture corrected, assuming a point source 
correction of 0.31 magnitudes. Here, z refers to zgso for the GOODS sources, 
and z' for the COSMOS sources. Non-detections are 2a. 



